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We calculate form factors of half-BPS operators in TV = 4 super Yang-Mills theory at 
tree level and one loop using novel applications of recursion relations and unitarity. 
In particular, we determine the expression of the one-loop form factors with two 
scalars and an arbitrary number of positive-helicity gluons. These quantities resemble 
closely the MHV scattering amplitudes, including holomorphicity of the tree-level 
form factor, and the expansion in terms of two-mass easy box functions of the one- 
loop result. Next, we compare our result for these form factors to the calculation of a 
particular periodic Wilson loop at one loop, finding agreement. This suggests a novel 
duality relating form factors to periodic Wilson loops. 
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1 Introduction 



One of the important challenges ahead is that of relating two apparently discon- 
nected realms, that of scattering amplitudes, and that of correlation functions. The 
first candidate theory to study is naturally maximally supersymmetric Yang-Mills 
theory (SYM). A step in this direction was taken in [THS], where correlation func- 
tions of composite operators in M = 4 SYM theory were considered in certain special 
lightlike limits where the distances between adjacent insertion points become null. 
Surprisingly, it was found that in the lightlike limit considered there is a simple re- 
lation between (n + /)-point correlators with I insertions of the Lagrangian, and the 
integrand of the n-point MHV amplitude evaluated at / loops. 

In this paper we will focus on an interesting class of physical observables which sit 
between the two worlds of amplitudes and correlation functions. These are the form 
factors, i.e. matrix elements of gauge-invariant, composite operators O between the 
vacuum and some external scattering stateQ 

F(l,...,n) = (l,...,n|O(0)|0> . (1.1) 

The external state can in principle contain all different particles in the theory. At 
strong coupling, these objects have been considered recently in [HE], where it was 
found that their calculation is technically equivalent to that of a periodic Wilson loop 
whose contour is specified by the lightlike momenta of the scattered particles. In 
distinction to the amplitude calculation at strong coupling [B], the momenta do not 
sum to zero as there is an operator insertion carrying momentum — q; furthermore, 
the sum of the momenta, ^2^ =1 Pi = q is not null, q 2 ^ 0, and defines the period of 
the Wilson loop contour relevant at strong coupling. It would be interesting to see 
whether this Wilson loop/form factor connection can be extended to weak coupling 
as was the case for the Wilson loop/amplitude duality [SHE]. In this paper we will 
find some evidence for this at one loop. 

At strong coupling, and leading order in 1/ \/A, the form factor calculation is 
insensitive to the polarisation of the external particles and, importantly, to the precise 
choice of the operator, as long as its anomalous dimension is small compared to Va [5] . 
This is similar to a feature of the amplitude calculation of [6], where the polarisations 
of the external states are only expected to play a role starting at one loop in the 
1 / vA expansion. At weak coupling, we have no reason a priori to expect that the 
form factor will be independent of the choice of the operator. We will then focus on 
the simplest class of operators in M = 4 SYM, that of half-BPS protected operators. 

To be specific, we will consider the operator Tr(0i20i2) inserted between a state 

1 We note that the dimension of the form factor in (jl.ll) is equal to d@ — n, where da is the 
physical dimension of O. 
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containing two scalars and an arbitrary number of gluons with positive helicity. The 
simplest member of this class is constructed with an external state containing just the 
two scalars, and is the Sudakov form factor F(q 2 ) : = (0i2(pi)0i2(p2)|Tr(0i 2 0i 2 ) (0) |0) , 
where q := pi + P2- This quantity was calculated at one and two loops in Af = 4 
SYM in a pioneering paper [5] by van Neerven, using an approach based on unitarity 
applied to Feynman diagrams developed in [TO]. Remarkably, the two- loop Sudakov 
form factor was found to satisfy an iterative relation very similar to that discovered 
in [II1EI2] for the four-point MHV scattering amplitudeso 

We begin in Section [2] by rederiving the one-loop result for the Sudakov form 
factor with an elementary application of unitarity. In Section [3] we will move on to 
consider a more generic class of n-point form factors containing n — 2 positive-helicity 
gluons in addition to the two scalars. We will first determine their expression at 
tree level for arbitrary n using recursion relations. Our result, Equation (13.31) . is a 
holomorphic function of the spinor variables describing the particle momenta, and is a 
close relative of the Parke- Taylor MHV scattering amplitude. For example, it localises 
on a complex line in twistor space, as was found in [13] for the MHV amplitude. 

We will then address the one-loop calculation of this class of form factors con- 
structed from a BPS operator and a state containing two scalars and an arbitrary 
number of gluons. Our strategy will consist in applying unitarity directly at the 
level of the form factor, thus bypassing Feynman diagrams, in the same spirit as the 
unitarity-based approach of [Hl[T5] for scattering amplitudes. The quantities enter- 
ing the form factor cuts are tree-level form factors and amplitudes. Both are gauge 
invariant objects, whose compact expressions we will recycle inside the loops, thereby 
obtaining simple one-loop integrands. The final result for this class of form factors at 
one loop is remarkably simple - see Equation (13.171) . It also very reminiscent of the 
expression for an n-point MHV amplitude at one loop, as we will describe. 

This close similarity will serve as an inspiration for Section HJ where we compare 
our weak-coupling result (13.171) for the form factors to a one-loop calculation of the 
periodic Wilson loop which, at strong coupling, computes the same quantity. Quite 
surprisingly, we will find that this periodic Wilson loop calculates the form factor 
also at weak coupling, or more precisely, the ratio of the one-loop form factor to its 
tree-level expression. Some issues related to gauge invariance of the prescription will 
also be discussed. 

Finally, we present our conclusions and discuss directions for further research in 
Section [51 A few appendices complete the paper. In particular, in Appendix [A] we 
write down a BCFW recursion relation [THldT] for the n-point tree-level form factor 
discussed earlier. In Appendix [B], we review for the reader's convenience the iterative 

2 We review this iterative structure of the form factor in Appendix 151 
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structure for the Sudakov form factor in Af = 4 SYM found in [9J. Appendix O 
contains the expression of the finite two-mass easy box functions, which appear in 
the form factor and Wilson loop calculations. 



2 The Sudakov form factor 



We begin by introducing the operator 

Oabcd ■= Tr(4> AB 4> GD ) - ^e ABCD Tr(0 LM LM ) , (2.1) 

where (j) AB := (l/2)e ABCD (j)cD- The operator in (12. II) belongs to the 20' representation 
of the SU(4) .R-symmetry group and is half BPS, i.e. it has vanishing anomalous 
dimension. Without loss of generality we will focus in the rest of this paper on its 
particular component 

O := Tr(0 12 12 ) . (2.2) 

In this section we concentrate on the simplest form factor one can construct using 
this operator, namely the two-point or Sudakov form factor, 

F{q 2 ) ■= (<Mpi)<Mp 2 )|0(O)|O) , (2-3) 

where 

q ■= P1+P2 ■ (2.4) 

As mentioned in the Introduction, this form factor in M = 4 SYM was studied at one 
and two loops in [9] , and calculated in [12] up to three loops. Here we will reproduce 
the one-loop result of [9] as a simple application of unitarity. We have also performed 
a similar calculation at two loops and found agreement with the result of [9], but we 
will not discuss it here. 



2.1 The one- loop Sudakov factor from unitarity cuts 

The Sudakov form factor is very simple to compute since, by Lorentz invariance, (12. 3p 
depends only on q 2 \ it can then be determined entirely from its unitarity cut in the 
q 2 channel. In [9j[l0] this procedure was carried out at one and two loops; here we 
depart from the approach of these two references in that we apply unitarity directly 
at the level of the form factor, bypassing Feynman diagrams. 
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(a) (b) 

Figure 1: In Figure (a) we show the diagram calculating the cut in the q 2 -channel of 
the Sudakov form factor (12. 3D . The cross denotes a form factor insertion. A second 
diagram with legs 1 and 2 swapped has to be added and doubles up the result of the first 
diagram. The result of this cut is given by (twice) a cut one-mass triangle function, 
depicted in Figure (b). 



The g 2 -cut of the form factor (i.e. its discontinuity in the g 2 -channel) is obtained 
from the diagram on the left-hand side of Figure HJ whose expression iqj 



F(1 V)| g2 _ cut = 2 y^LIPS(/ 1 ,/ 2 ;g)F( o )(/ 1 ,/ 2 ;g)A(°)(0 12 (p 1 ),0i2(p2),034(/i),034(/2)) 

(2.5) 

where the Lorentz invariant phase space measure is 



dLIPS(Zi, l 2 ; q) := d u h d u l 2 5 + (li)5 + (l 2 2 )5 u (h + l 2 + q) , (2.6) 

and q is given in (|2.4p . The tree- level component amplitude appearing in (12. 5p . 
(4>i 2 (pi), <Pi 2 (p 2 ), ^34(^1), ^34(^2)) , can be extracted from Nair's superamplitude 



Amy := 9 n ~ 2 (2n)^ ( £ A,Aj ( £ A^) f[ — , (2.7) 

i=l i=l i=l \ U 1 

where A n+ i = Ai. The result is 

A^^l 2 {pi)Al2{P2)Au{h\ Mk)) = ■ (2.8) 

The other quantity appearing in (1231) . is the tree-level expression for the form 
factor (12.31) . which is trivially equal to 1. Thus, we get 

(12) (hh) _ „ 2 f t-i- TT-)ri / 7 , v 1 



F(1) (g 2 )| g ^ cut = 2 JdLIPS(h,l 2 ;q) ™^ = -2q 2 JdLIPS(k,l 2 ;q) 



(h+Pi) 2 ' 
(2.9) 



3 In this and the following formulae we omit a power of the 't Hooft coupling, defined as a :- 
(3> 2 7V)/(167r 2 )(47re- 7 ) e . Note that this is 1/2 the 't Hooft coupling defined in QJ]. 
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The expression in (12.91) represents the cut of a one- mass triangle T lm (q 2 , e), depicted 
in Figure [H(b). It can instantly be lifted to a full loop integral since it depends only 
on a single kinematic invariant, much in the same spirit as [19] . Doing so we get 

F«(g 2 ) = 2Tri(g 2 ,e) , (2.10) 



where the function Tri(g 2 ,e) := q 2 T lm (q 2 ,e) is explicitly evaluated in (1B.3j) . This 
result agrees with that of [9]. 



3 Multi-point form factors 



With the operator O introduced in (12.21) . we can now construct an infinite sequence 
of n-point form factors, 

F(l, . ... ,n;q) := (g + (pi) ■ ■ ■ MPi) ■ ■ • <Mpj) ■ ■ ■ g + (p n )\0 (0)\0) , 

(3.1) 

i.e. we take matrix elements of O between the vacuum and a state containing the 
same two scalars already appearing in (I2.3p . along with n — 2 positive-helicity gluons. 
The particular form factor we consider is colour-ordered with respect to the positions 
of the external particles. However notice that the operator O is a colour singlet, 
hence the momentum q it carries can be inserted at any position in colour ordering. 
We now present the calculation of (13.11) . first at tree level and then at one loop. 



3.1 Tree level 



At tree level, it is easy to calculate the form factor (13.11) . Indeed, we observe that 
factorisation theorems are valid also for form factors - in fact they apply to Green's 
function in general, see for example [20J for a discussion. We can then use tree- 
level factorisation in order to write down a BCFW recursion relation [T6l[T7] for form 
factors. Our result for this quantity is very simpleQ 



^(1 ife,-,^ = ^- 2 (2vr)V 4 )(^A fc A fe -g) F MHV , (3.2) 

where 



k=l 



MHV 



(12) ••>!> 



(3.3) 



4 The calculation of (|3.2j) is presented in Appendix [5] We have also checked our result against 
Feynman diagrams in a few cases. 
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Here p m := X m X m , and YlZi=iPrn '■= Q is the momentum carried by the operator 
insertion. A number of remarks are in order. 

1. The expression in (13. 3 P is purely holomorphic in the spinor variables of the 
external particles. In this sense this form factor is the closest off-shell relative of the 
Parke-Taylor MHV amplitude. We will refer to this form factor as to the "MHV form 
factor". Note however that the momenta do not sum to zero. We also note that 
this expression is very reminiscent of the formula for the infinite sequence of Higgs + 
gluons amplitudes considered in [2Tj . 

2. We can easily transform the tree-level form factor (13 .2\i to Penrose's twistor 
space, with the result 



/n 
d 4 x \\5^^ m + x\ 

m=l 



(3.4) 



where (A m , /i m ) are the twistor space coordinates of the m th particle. As for the 
case of the MHV amplitude, holomorphicity of the tree-level form factor ( 13 .3p ensures 
that -Fmhv can be pulled out of the half-Fourier transform. Equation H3.4[) shows that 
our MHV form factor is localised on a line in twistor space, similarly to the MHV 
amplitude 



3. We notice that (13. 3p satisfies an auxiliary (helicity) condition 

o(- X m ^ + ~ X m-^)F(l,...,n) = h m F(l,...,n) , (3.5) 

where h m is the helicity of the m th particle. This relation is just the statement that 
each external state must transform appropriately under the little group of a massless 
vector [TBI. 



3.2 The one-loop MHV form factor from unitarity 

We can now insert the compact expressions for the tree-level form factor (13. 3p found 
earlier into unitarity cuts, and glue it with tree-level amplitudes in order to build 
cuts of loop form factors. Our strategy will consist in calculating the cuts of form 
factors in all kinematic channels, and then reconstructing the function which has all 
the correct cuts, in complete analogy with the approach of [H1[T3] for amplitudes. 
The final result will be expressed in terms of the elements of the basis of one-loop 
scalar integral functions, times rational coefficients!! 

5 This approach might of course miss purely rational terms, not linked to discontinuities of integral 
functions. Since we work in M = 4 super Yang-Mills such terms are absent. 
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The operator we have chosen in (12. ip is unrenormalised; as a consequence, the 
form factor will contain no ultraviolet-divergent functions. At one loop this implies the 
absence of bubbles. Unlike amplitudes, we will encounter one- mass triangles; however, 
as we shall see momentarily, their purpose will be that of canceling certain otherwise 
unwanted infrared divergences in multi-particle channels. Indeed, anticipating our 
story a little, the result of our calculation will consist of a sum of infrared- divergent 
terms containing only two-particle kinematic invariants, along with a sum over all 
possible finite parts of two-mass easy box functions. 

We now discuss all possible cuts of the form factor. Here is an outline of the main 
features of our calculation: 

1. We will perform cuts in three distinct channels. The first one is the q 2 - 
channel, where the tree-level form factor entering the cut is a Sudakov form factor. 
In the second case, which needs to be treated separately, a tree form factor with one 
additional external particle enters the cut. These two cases give rise to two- mass 
easy box functions, as well as to triangles. Finally, in the most generic channel, the 
form factor entering the cut expression has an arbitrary number of legs. Importantly, 
in this case we will discover that, after factoring out the tree-level expression of the 
form factor, the integrand will be the same as that of the one-loop MHV amplitude 
in M = 4 super Yang-Mills, calculated in [Hj. 

2. In principle we will have to perform sums over all possible internal states which 
can run in the loop. We have found that, after performing the sums over different 
states when necessary, the cut-integrand becomes the same in all channels. This 
is in complete analogy with the calculation of [T4] of the one-loop MHV scattering 
amplitude in M = 4 SYM. 

3. The form factor is ordered with respect to the position of the external particles 
in colour space. One may be tempted to think of the form factor as an amplitude 
with one (or more, if there are more operator insertions) leg going off shell, but this 
picture would not be true as far as the colour ordering is concerned since the operator 
insertion is a colour singlet. Hence, its position - the position of the "off-shell line" - 
does not affect the colour ordering of the external states. One must therefore insert 
the operator in all possible ways and sum over the corresponding contributions. In 
practice, this possibility will arise only in the g 2 -cut diagrams. 
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3.2.1 The q 2 channel 



We begin by computing the cut in the q 2 channel, where 

n 

q = , (3.6) 
k=i 

is the sum of the momenta of the external particles. This cut is represented in Figure 
El and is given explicitly by 




Figure 2: The q 2 -cut of the one-loop form factor. Note that the complete cut is 
obtained by summing over cyclic permutations of (1, ... ,n). 

FW(l,...,n; g )| g2 _ cut = fdUPS(l u h;q) £ (3.7) 

J V{l,...,n) 
[F^(h, Z 2 ; q) A® (<f>M, hM,g + (pi), MPi), • • • MPj), g + (p n )) } ■ 

Note that we have included a sum over cyclic permutations of the particles 1,2, ... ,n. 
This is because the insertion of the operator does not affect the colour ordering, 
therefore all terms in the sum have the same ordering in colour space. 



We now move on to describe the ingredients appearing in (13. 7p . The tree-level 
form factor is equal to 1. The tree amplitude A^°\ obtained from (12. 7p . is 

(Z 2 i)(i2)---K)(y 2 ) • 1 • ] 

Let us focus on one term in the sum, specifically that with the ordering l\, l 2 , 1, . . . , n 
of the legs. The others will be obtained by cyclically permuting 1, . . . , n. This term 
can be written as 



F(°)(l,...,n; q) JdLTPS{h,l 2 ;q) 



{kl){nh) 
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where F^°'(l, . . . , n; q) is the tree-level MHV form factor ( 13.31) . It is easy to show that 



(W(nl) 
(hi) (nix) 



2 (piq)(hp n ) ~ (pih)(qp n ) + (piPn)(qh) 



(k+Pl) 2 (ll+Pn) 2 

2(qPn)(qpi) - q 2 (piPn) . (piq) (qpn) 



(3.9) 



(h+Pl) 2 (h+Pn) 2 (h+Pl) 2 (h+Pn) 



2 ' 



where (p a Pb) '■= Pa'Pb- This reduction thus leads to the g 2 -cut integral of a two- 
mass easy box and two scalar triangles, all of which can be calculated using standard 



box obtained here has massless legs p\ and p n , whereas the momenta of the massive 
corners are — q and P^n-i = 1 ~ Pn — Pi, where 



Furthermore, we have two two-mass triangles. The first one has massless leg p\ and 
massive legs — q and q — Pi] for the second we just replace p\ by p n . Notice that 
by summing over all permutations we will generate twice every possible two-mass 
triangle, where the massless leg is in turn any one of the momenta of the particles, 
and one of the massive corners has momentum —q. Furthermore, we will produce all 
possible two-mass easy boxes with massless legs a — 1 and a, and massive legs — q and 
Pa+ia-i for any a € (1, . . . ,n). 

Rather than calculating these dispersion integrals, we will now proceed to inspect 
other kinematic channels. This will enable us to reconstruct the one-loop form factor 
from its cuts. 



3.2.2 The cut in a generic kinematic channel 



Figure 3: The cut of the one-loop form factor in the s a +i.&-i- channel. 

Here we discuss the most generic cut, depicted in Figure [3) There are three possibili- 
ties to consider. The two scalars can either be emitted both from the form factor; or 

6 The important observation that dimensionally regularised dispersion integrals are well-defined, 
and rather simple objects to compute, was made in |10j . 




Pab '■= Pa H VPb ■ 



(3.10) 
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both from the MHV amplitude; finally, one can be emitted from the form factor, and 
one from the amplitude. In the latter case, it is necessary to sum over gluons and 
fermions, which can propagate in the loop. Remarkably, the three cases give rise to 
the same integrand. For this reason, in Figure |3] we need not distinguish the different 
types of particles. 

For the sake of definiteness, we focus on the case where two scalars are emitted 
from the amplitude. The corresponding expression is quickly seen to be 

F { \ a+ib _^ ut = JdLlPS(h,l 2 ;P a+lb ^) (3.11) 

F(°\-li,-li,V + ,...,a* + ; q ) 

A^(lf,lt,(a + iy + ./" (!>■ I)"') , 

where P ab is defined in (13.101) . and the superscript indicates the particle's species. 

Notice that the particles emitted from the tree-level form factor fix the colour 
ordering, hence there is no sum over permutations to be performed in this case. 
Using the explicit expressions of the tree-level form factor and the amplitude, (13. lip 
becomes 



(aa + l)(kh) (b-lb)(kl 2 ) 



{al 2 )(l 2 a + l) (b-lh){kb) ' 

(3-12) 

where is now the tree- level form factor (13. 3p . The reader may have recognised 
that the integrand of (13. 12[) is the same as that appearing in the two-particle cuts of 
the one- loop MHV amplitude considered in [H] . Performing the reduction as in that 
paper, one finds that the integrand of (13.121) can be recast as the sum of four terms, 

R(b, a + l) + R(b- 1, a) - R(b, a) - R{b - 1, a + 1) , (3.13) 

where 

(bl 2 )(ali) 2 (hPb)(hPa) + (hPa)(hPb) - (hh)(p a Pb) 

R(b, a) : 



(bh)(al 2 ) 2(hPb)2(hp a ) 

1 (PbP) (PaP) 2(p a P)( Pb P)-P 2 (p aPb ) 

+ 2( Pb h) + 2(pJ 2 ) + 2(p b h)2(pJ 2 ) ' 1 • ] 

and we used momentum conservation l\ + 1 2 + P — 0, where P is the cut momentum 
(P 2 = s a+ \ ib -i in this cut). We also used the cut conditions l\ = l\ = 0. 

Equation (I3.14p shows that each i?-function gives rise to a bubble, two triangles 
and a two- mass easy box function. When all four of the terms in ( 13. 13[) are present, 
bubbles and triangles cancel. We will see in the next section that when b = a there are 
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Figure 4: On the left, we represent a two-mass easy box function. The momenta 
p a and pi, are null, whereas, in general the remaining momenta P := P a+ ib-i and 
Q := — q + P(, + i a _i are not null. The cases when either P 2 or Q 2 , or both, are 
also null, correspond to the one-mass and zero-mass boxes, obtained as smooth limits 
from the expression fIC.ll) of the two-mass box function. On the right, we represent 
a two-mass triangle arising from the cuts considered in Section l3.2.3[ The thick line 
represents the momentum carried by the operator. 



some surviving triangle functions (whereas bubble always cancel among themselves, 
as anticipated). 

Specifically, the box function arising from R(b, a) has massless legs a and b, and 
massive legs P a +ib-i, and — q + Pb+\ a -i- Furthermore, it appears here in the corner 
cut s a+ i^-i, i.e. the cut momentum legs are those adjacent to the corner momentum 
-Pa+16-i- The other .R-terms in ( 13 . 13f) will give the same Sa+i^-i-cut of different two- 
mass easy box functions, where the entries I and m of R(l, m) denote the massless legs 
of the box. Note that the same box in (13.1 4ft will appear in all of its other cuts, with 
precisely the same coefficient. To see this, notice that the quantity 2(p a P){pi > P) — 
P 2 {j>aPb) is invariant under P — > P + ap a + (3pb for any a and (3. We remind the 
reader of the pictorial representations of the box and triangle functions in Figure HI 



3.2.3 The cut in a (q — p a ) 2 channel 



Next we consider a channel where a single particle, of momentum p aj is emitted from 
the form factor - see Figure [51 This (q — p a ) 2 -channel cut can be regarded as the 
special a = b case of the generic s a+ i ;,_i-channel cut described in the previous section. 
Note that in the corresponding calculation of the one-loop MHV amplitude, with a 
single particle of momentum p a emitted from a vertex, this cut would be absent since 
p 2 = 0. However, the momentum flowing from the form factor is here — q + p a , which 
is not null, and there is a new nontrivial cut to compute. 
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Figure 5: The (q — p a ) 2 -cut of the one-loop form factor. 



The cut integral is given by a formula identical to f)3.12p except that now a = b. 
The result for the cut is therefore 



R(a, a + 1) + R(a — 1, a) — R(a, a) — R(a — 1, a + 1) 



(3.15) 



where, using (I3.14p . we see that R(a,a) = 1. As a consequence, bubbles still cancel 
among the four terms, but there are uncanceled triangles (as well as boxes). Specif- 
ically, using the second line of f l3.14p . it is immediate to see that there is a single 
surviving triangle, which gives a contribution 



(PaP) 
2(Pah) 



o+ln-1- 



(3.16) 



This is therefore a two- 



where in this particular cut, P = —q + p a = 
mass triangle, where the massless leg has momentum p a , and the two massive legs 
have momentum P a +i a -i and — q, respectively. We have thus found the same set of 
two-mass triangles as we found from the g 2 -cuts discussed earlier. 



3.3 The complete result 



The next step consists in inspecting the integral functions appearing in all cuts; 
one then reconstructs the function which has the correct cuts in all channels. An 
important observation is that the triangle functions (with the precise coefficient they 
appear with) precisely cancel all infrared-divergent terms in the q 2 and (q — p a ) 2 
multi-particle channels. The remaining multi-particle infrared-divergent terms cancel 
among all the boxes in the same way as in the calculation of the one-loop MHV 
amplitude. 

After these cancellations are performed, the final result (divided by the tree-level 
form factor) has a remarkably simple form: it consists of the sum of infrared- divergent 
terms containing only two-particle invariants made of momenta which are adjacent 
in colour space, plus a sum of finite parts of two-mass easy box functions, where the 
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massless legs are any two of the particle's momenta, all with coefficients equal to onelll 
F^(l,...,n ]q ) = F(°\l,...,n;q)[-J2 ,2 + Fin 2m > a ,p b , P, Q) • 

1=1 a,b 

(3.17) 

In the sum on the right-hand side of ( 13. 17ft the external particles are distributed fol- 
lowing cyclic ordering. Note however that the operator insertion, carrying momentum 
—q, does not affect this ordering. 

As an example, we list below the finite box functions appearing in the four-point 
result: (1, {-q}, 2, {3, 4}) (massless legs 1 and 2); (1, 2, 3, {4, -q}), (1, {2, -q}, 3, 4) 
(massless legs 1 and 3); (1, {2, 3}, 4, {— q}) (massless legs 1 and 4); (2, {— q}, 3, {4, 1}) 
(massless legs 2 and 3); (2, 3, 4, {1, — q}), (2, {3, — q}, 4, 1) (massless legs 2 and 4); 
(3, {— q}, 4, {1, 2}) (massless legs 3 and 4). Here we denote by (a, {P}, b, {Q}) a two- 
mass easy box with massless momenta p a and pb, and corner momenta P and Q. 



4 Comparing form factors to periodic Wilson loops 



In [Ul5] , a prescription to calculate form factors at strong coupling using the AdS / CFT 
correspondence was proposed. We recall that the form factor calculation, as well as the 
scattering amplitude calculation at strong coupling [4ll6 pTH30] , are both equivalent to 
the problem of computing minimal surfaces in AdS space, and hence Wilson loops, but 
the boundary conditions in the two cases are quite different. For n-point amplitudes, 
the boundary (the contour of the Wilson loop) is the n-edged closed polygon obtained 
by joining the lightlike momenta of the particles following the order induced by the 
colour structure of the planar amplitude. For n-point form factors, the boundary is 
an infinite periodic sequence of n lightlike segments [HIS], see Figure E] for an example 
of a form factor with three particles, of momenta p\, P2 and p^. The period q is the 
momentum of the inserted operator, which at strong coupling corresponds to having 
an additional closed string state inserted on the worldsheet. Therefore, in the form 
factor case the worldsheet stretches all the way from the boundary at r = to r — > oo, 
where r is the radial coordinate of the T-dual AdS space (HE]- 

It was observed in [7JIE] that the one-loop calculation of a Wilson loop with the 
same polygonal contour as in the strong-coupling calculation reproduces the one-loop 
MHV amplitude, divided by the tree amplitude. This was the first manifestation of a 

7 A more explicit way to write the sum of finite two-mass easy box functions in ([3.1 7p is 
EILi Ella -1 (i - (l/2)£„/ 2 -i,r) Fin^°, where the relation to the functions Pin 2me (p a ,Pb,P,Q) 
depicted in Figure[4]is obtained by setting p a = Pi-i, pt — Pi+r, and P = pi + ■ ■ ■ +Pi+ r -i- We give 
the explicit expression of these finite box functions in Appendix O 
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9 



Figure 6: Periodic Wilson loop featuring in the duality with form factors. In this case 
the period is q = p\ + p 2 + P3 ■ 



new amplitude/ Wilson loop duality at weak coupling, which was further checked and 
studied in [3TH4"2] . An immediate question is therefore whether there is also a form 
factor/periodic Wilson loop duality at weak coupling, i.e. whether the form factor 
can be equivalently evaluated at weak coupling from a Wilson loop with the same 
periodic contour found at strong coupling. In this section, we calculate the periodic 
Wilson loop at one loop, and compare it to the result (13. 17ft we obtained earlier for 
the n-point form factor using unitarity. We will find that the two are in agreement. 



The Wilson loop we consider is 



W[C n 



Tr V exp 



ig / dr x^(t)A^(x(t)) 



(4.1) 



where an example of a periodic contour C n with n = 3 is depicted in Figure El For the 
calculation it is convenient to choose the Feynman gauge and, as discussed in [7J|8], 
this leads at one loop to two classes of diagrams, namely the infrared- divergent cusp 
diagrams, where a propagator connects two adjacent edges, and those diagrams where 
the two edges connected by a propagator are not adjacent. Diagrams in the latter 
class are finite. 



The first evidence of the duality is straightforward to detect, and comes from the 
infrared-divergent part. The form factor result (13.171) shows that the only divergent 
terms contain two-particle kinematic invariants made of momenta that are adjacent 
in colour ordering, 

F^(l,...,n;q)\ IR = F<°>(1, . . . , n; q) 

with s ii+ i := (pi + Pi+i) 2 - This contribution exactly matches the sum of the one- 
loop cusp diagrams for the Wilson loop, in an identical way as for the matching of 
the infrared- divergent parts of one- loop MHV amplitudes and Wilson loop 0[8]|^| 
As in that case, the dual Wilson loop calculation does not reproduce the tree-level 

8 We recall that a cusp contributes a term — (— s ii+ i) _e /e 2 [71151 



1 - 

'72 



~ s ii+lj 



(4.2) 
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Figure 7: We illustrate in this figure the correspondence between Wilson loop diagrams 
connecting non-adjacent edges and finite parts of two-mass easy box functions. Notice 
that this correspondence holds in the Feynman gauge. 




(a) (6) (c) 



Figure 8: One-loop finite Wilson loop diagrams in the three-point case. Diagrams 
(a), (b), (c) are equal to the finite parts of the two-mass easy boxes (1, 2, 3, {— q}), 
(2, 3, 1, {—?}), and (2, {— q}, 3, 1), respectively. 

prefactor. Note that each period contains n cusps, and we sum once over these n 
cusps. In particular, for the two-point case, the period contains two cusps which give 
an identical contribution. This is related to the factor of two in the Sudakov form 
factor gHEE] in fl2TT0|) . 

Next, we move to the finite parts. We found in ( 13.1 7p that the one-loop form 
factor contains a sum of finite parts of two-mass easy box functions (and one-mass 
boxes, as special cases of the former). It was shown in [8] that a Wilson loop diagram 
where a propagator connects two non- adjacent edges, called p and r on the left-hand 
side of Figure [71 is identical to the finite part of a two-mass easy box function with 
massless legs p and r, right-hand side of the same figure. The massive corners P 
and Q of the two- mass easy box are mapped to sums of adjacent momenta on the 
Wilson loop side, as shown in Figure [71 This crucial observation puts in one-to- 
one correspondence the finite boxes appearing in (I3.17P with all the one-loop Wilson 
loop diagrams connecting non-adjacent edges in one period. For example, the three- 
point form factor contains the finite parts of three two-mass-easy boxes (1, 2, 3, {—q}), 
(2, 3, 1, {— q}), and (2, {— q}, 3, 1). These functions are exactly given by the three 
Wilson loop diagrams depicted in Figure [SJ It is not difficult to see that this mapping 
is true for general n-point form factors. 

A few remarks are in order here. 
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Figure 9: In this figure we represent a one-loop Wilson loop diagram that should not 
be included. 

1. Notice that q is not part of the contour, and is therefore not connected to any 
other edge by propagators. 

2. We draw all finite diagrams where we connect all pairs of non-adjacent legs 
within a period. We can in all cases map these configurations to finite two-mass easy 
box functions using the correspondence of Figure [7J 

3. In principle one may think that the full calculation of a periodic Wilson loop 
would require a summation over all translationally inequivalent diagrams. In partic- 
ular one should also include diagrams where a propagator stretches over more than 
one period, as for example that in Figure |9j As we have seen, in order to reproduce 
the form factor result these diagrams should not be included. This means that we 
need to make a truncation of the periodic Wilson loop and only consider the diagrams 
within one period. Notice that at strong coupling the form factor is calculated by the 
area of one period [HE]. 

4. By restricting to one period, we have succeeded in mapping the truncated 
periodic Wilson loop calculation, which is not obviously gauge invariant, to the form 
factor, which is gauge invariant by definition. It would be important to explore this 
issue further, in particular through higher-loop calculations of the periodic Wilson 
loop. We leave this for future work. 

5. An important difference between the amplitude and the form factor calculations 
at strong coupling is the presence in the latter case of additional boundary conditions 
at r — > oo due to the operator insertion on the worldsheet. It would be important 
to understand what modifications of the weak-coupling Wilson loop calculation this 
entails. This might require the insertion of additional operators compensating the 
gauge non-invariance of the calculation^ 

9 We thank Gregory Korchemsky for discussions on the issue of gauge invariance. 
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5 Conclusions 



To summarise, we have calculated form factors of half-BPS operators in M = 4 super 
Yang-Mills theory at tree level and one loop. This took advantage of the powerful 
methods of recursion relations p^6|[T7] and unitarity [i~4"l[l~5] . which may also be applied 
to form factors. The expressions we found for the one-loop form factors with two 
scalars and an arbitrary number of positive-helicity gluons in particular share salient 
features with MHV scattering amplitudes. We then found that these form factors 
could be straightforwardly derived from periodic Wilson loops, suggesting that there 
is a new duality between form factors and periodic Wilson loops. 

One can also easily generalise the form factor defined in (13. ip to one containing 
external states with different particles, in particular negative-helicity gluons. These 
should correspond to non-MHV extensions of the form factor we considered. It would 
be useful to explore a supersymmetric formulation of such quantities. Most impor- 
tantly, it would be very interesting to consider form factors with insertions of different 
operators, including non-BPS ones such as the Konishi operator. 

We also found that the tree-level form factors studied here localise in twistor 
space, and a similar question arises as to whether more complicated form factors 
(e.g. with insertions of the Konishi operator) enjoy similar localisation properties 
in twistor space. We anticipate that form factors similar to (13. ip . but where the 
external states also contain negative-helicity gluons, will be localised on unions of 
lines as found in [T3] for the tree-level non-MHV scattering amplitudes. One would 
also like to investigate the application of the whole MHV diagram approach to form 
factors, much in the same spirit as was done in [21] when considering interactions 
with massive bosons. 

Finally, an obvious important question is to investigate this conjectured form fac- 
tor/Wilson loop duality at higher loops. It seems likely that there will be many 
further interesting and productive directions to explore as part of applying the tech- 
niques and insights obtained from the study of scattering amplitudes to correlation 
functions more generally. 
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A Recursion relations for tree-level form factors 



In this appendix, we apply the BCFW recursion relation to calculate the tree-level 
form factor defined in (12. 3p . The basic inputs are three-point amplitudes, plus the 
three-point MHV and anti-MHV form factors, 

F 3 iH^M- {jk){ki) > F s ^^ k ^-JjJ^-y ( A - x ) 
which can be easily derived by using Feynman diagrams. 

As is standard in the BCFW recursion relation, we use an shift, Aj — » Aj — z\j, 
Xj —> \j + z\i. Thus we obtain a one complex parameter family of form factors, 
F(z). As mentioned earlier, factorisation theorems are also valid for form factors; 
therefore, exactly as in the case of scattering amplitudes, by using Cauchy's theorem 
we can calculate the form factor by summing the residue of the poles from various 
factorisation channels (we also require F(z) — > as z — > oo), 




as shown in Figure [TDJ Notice that there are two diagrams, since we can insert the 
operator either on the left- or on the right-hand side. 




Figure 10: BCFW recursion diagrams for tree-level form factors. 

As an example, we will derive the expression for the MHV form factor. Since 
we require the form factor F(z) to vanish as z — > oo, we shift the momenta of two 
positive-helicity gluons. It is not difficult to see that the non-vanishing diagrams 
are those where the right-hand side is an anti-MHV three-point amplitude. Without 
loss of generality, we calculate the n-point form factor F n (l^, 2 + , 3^, 4 + , ■ ■ • , n + ) by 
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performing a [2,4) shift. There are two non-vanishing channels to consider. The 
contribution of the first one is 



F n -i(5 + , ■ ■ ■ ,n + , l^, 2 + , P^ 12 ) — A^P^, 30 12 ,4 + ) 

_ (1 P) 2 1 [3 4][4P] 

" (5 6) •••(2 P)(P 5)s 3 4 [P 3] ' 

where P = P 34 + Z34A2A4, z 34 = —(4 3)/(2 3). The other one is 

F n _!(6 + , • • • , i> 12 , 2 + , 3 012 , P l+ )— A 3 (P'-, 4 + , 5+) 

S45 

_ (1 3) 2 1 [4 5] 3 

" (6 7) • • • (3 P')(P' 6) s 45 [P' 4] [5 P'] ' 



(A.3) 



(A.4) 



where P' = P45 + Z45A2A4, Z45 = —(4 5)/(2 5). The sum of these two terms indeed 
reproduces the MHV form factor result (I3.3p . Notice that as an illustration we chose 
a complicated shift; it would be much simpler to shift the momenta of two adjacent 
gluons. 



B The iterative structure of [9] for the two-loop 
form factor 



In [9], Van Neerven calculated the one- and two-loop contribution to the Sudakov 
form factor F(q 2 ,e), defined in (13. ip . His results ari 



oop 
& 



F {l \q\e) 



2Tri(g 2 ,e) , 
4LT(g 2 ,e) + CT(g 2 ,e) 



(B.l) 
(B.2) 



where the one- loop one-mass triangle Tri(g 2 , e) and the two- loop ladder and crossed 
triangle, LT(g 2 ,e), CT(g 2 ,e) respectively, are given byf"1 



Tri(g 2 ,e) 



^_^ 7 r(l + e)r 2 (-e) 



r(l-2e) 



(B.3) 



2\-e 



e 2 2 3 s 1440 V 1 



10 In this appendix we indicate explicitly the dependence of the form factor on the dimensional 
regularisation parameter e. 

11 In the following formulae we actually divide these functions by a power of q 2 per loop. 
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LT(g 2 ,e) = {~q 2 )- 2e e 2 ^ 



-G(2,2)G 3 (2 + e,i,l) 



(B.4) 



-G(2,l) -G 3 (2,l,l + e)+G 3 (l,l,l) 
Le 



CT(g 2 ,e) 



,2\-2e 



„2\-2<e 



1 5tt 2 29 a 3 4 ^, . 
4^ + 24^ + 67 C3 + 32 7r + ° (e) 
1 7T 2 83 , 59 . 

? " F " 37 & " IS 1 + 0(e) 



where 



G(x,y) 
G 3 (x,y,z) 



T(x + y + e- 2)r(2 - e - x)T(2 - e - y) 

T{x)r(y)r(4 - x - y -2e) ' 
F(2-x-z- e)r(2 -y-z- e)T(-2 + x + y + z + e) 
T{x)T{y)T{A - x - y - z - 2e) 



(B.5) 

(B.6) 
(B.7) 



In [10], the functions LT(g 2 , e) and CT(g 2 , e) were derived from dimensionally-regulated 
dispersion integrals; here we present these functions in the form given in [52] . 

In [9] it was proved that the Sudakov form factor exponentiates at two loops. In 
formulae, one finds 



F {2 \q\e)- l -(F^\ q 2 ,e)) 



2\-2e 



% + ^ + 0(e) 



We can recast (1B.8P in an ABDK/BDS form, namely 
F^\ q \e)- l -{F^\ q \e)) 2 

where f$ (e) = f + /ie + f 2 e 2 . We then find 

/o = -2C 2 , A = -2Cs • 

(2) 

We also find a condition relating / 2 and Cp F , namely 



/ F 2 F ) (e)F( 1 )(g 2 ,2e) + 4 2 I ) + 0(e 



(2) _ h 

FF _ 2 IS 



(B.8) 



(B.9) 



(B.10) 



(B.11) 



On the other hand, the four-point MHV amplitude (divided by the tree-level ampli- 
tude) satisfies [TT|[T2] 



(B.12) 



M ( 2 )( e )-i(M«( e )) 2 = f^ 2 \t)M^\2t)+C^ + 0{t) 
with f^(e) = f + e f 1 + e 2 f 2 and 



fo = -Ca 



/o 



/i = -Cs 



/i 



(B.13) 
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The factor of 1/2 in the result ( IB. 131) is a matter of convention - it can be understood 
once one recalls that /o and f\ are written in a convention where the 't Hooft coupling 
obds is twice as that used in the present paper (as well as in [9]). Indeed, inspecting 
(IB. 121) one quickly realises that the combination af^ must be independent of any 
conventions used to define the coupling, since the left-hand side a 2 M^ is clearly 
convention independent. 



C The two-mass easy box function 



A compact form of the finite part of a two-mass easy box function containing only 
four dilogarithms was first derived in [53J, and then found independently in [22] in 
the context of MHV diagrams, where an analytic proof of its equivalence with the 
conventional expression of e.g. [54J was given. Expressing the two-mass easy box as 
a function of the kinematic invariants s := (P + p) 2 , t := (P + q) 2 and P 2 , Q 2 , with 
p + q + P + Q = 0, its finite part is 

Fin 2me (s, t,P 2 ,Q 2 ) = Li 2 (l-aP 2 ) + Li 2 (l-aQ 2 ) -Li 2 (l- as) -Li 2 (l-at) , (C.l) 

where 

P 2 + Q 2 -s-t 

a = P 2 Q 2 - st ■ (C - 2) 
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